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ABSTRACT 


An  analytical  study  was  undertaken  to  determine  two- 
dimensional  wall  conduction  effects  during  film  condensation 
on  an  internally  finned,  rotating  heat  pipe.  An  earlier 
Finite  Element  method  was  modified  to  solve  the  problem 
by  using  different  numbers  of  axial  increments  and  different 
numbers  of  elements.  Resulting  heat  transfer  rates  were 
shown  to  be  similar  to  those  obtained  earlier. 

A rotating  heat  pipe  was  tested  using  various  copper 
condensers  including  two  smooth-wall  cylinders  and  a 0.5 
degree  truncated  cone.  All  condensers  were  tested  at 
different  rotational  speeds  using  distilled  water  as  the 
working  fluid.  In  the  case  of  the  1.46  inch  diameter 
cylindrical  condenser,  data  were  also  taken  with  Freon  113 
and  ethanol.  The  heat  transfer  rate  of  each  condenser  was 
plotted  against  the  saturation  temperature  of  the  vapor. 

Using  the  experimental  results  obtained  in  this  thesis, 
together  with  the  data  of  earlier  theses,  a comparison 
was  made  with  the  theoretical  analysis  of  Roetzel.  Agree- 
ment is  reasonable,  and  Roetzel 's  analysis  can  be  used  to 
predict  overall  rotating  heat  pipe  performance  to  within 
± 20  percent  for  water.  The  data  for  Freon  113  and  ethanol, 
however,  fall  approximately  35  percent  higher  than  this 
theory . 
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I. 


INTRODUCTION 


A.  THE  ROTATING  HEAT  PIPE 

The  rotating  heat  pipe  is  a closed  rotating  container 
capable  of  transferring  large  amounts  of  heat  from  rotating 
machinery.  It  consists  of  three  main  parts:  (1)  a cylin- 

drical evaporator,  (2)  a cylindrical  or  truncated-cone 
condenser  and  (3)  a working  fluid.  The  operation  of  a 
rotating  heat  pipe  is  shown  schematically  in  Figure  1. 

When  the  rotating  heat  pipe  is  rotated  above  a critical 
speed,  the  working  fluid  forms  an  annulus  in  the  evaporator. 
Heat  input  to  the  evaporator  causes  some  of  the  working 
fluid  to  vaporize,  and  this  vapor  flows  into  the  condenser 
transporting  latent  heat  of  vaporization.  External  cooling 
of  the  condenser  allows  the  vapor  to  condense  on  the  inner 
condenser  wall.  The  centrifugal  force  due  to  the  rotation 
has  a component  acting  along  the  condenser  wall  which 
accelerates  the  condensate  back  to  the  evaporator. 

B.  ROTATING  HEAT  PIPE  LIMITATIONS 

Ballback  [1]  has  studied  the  limitations  imposed  on 
the  rotating  heat  pipe  by  various  fluid  dynamic  mechanisms 
in  a preliminary  theoretical  analysis  of  its  operation. 

By  using  existing  theories  and  experimental  correlations, 
he  was  able  to  estimate  the  sonic  limit,  the  boiling  limit, 
the  entrainment  limit  and  the  condensing  limit  of  performance. 
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Figure  1.  Schematic  Drawing  of  a Rotating  Heat  Pipe 
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A typical  rotating  heat  pipe  with  the  physical  dimen- 
sions as  shown  in  Table  1 was  chosen  in  order  to  demon- 
strate the  predicted  performance. 

Table  1 

Specifications  of  a Typical  Rotating  Heat  Pipe 


Length 

14.000 

inches 

Minimum  diameter 

2.000 

inches 

Wall  thickness 

0.125 

inches 

Internal  half  angle 

1.  000 

degree 

Rotating  speed 

2700 

RPM 

1.  Sonic  Limit 

As  the  heat  flux  in  the  rotating  heat  pipe  is 
increased,  it  is  possible  that  the  flow  rate  of  the  vapor 
down  the  container  may  be  limited  by  the  attainment  of  a 
choked  flow  condition  for  the  vapor.  This  choking  of  the 
vapor  flow  will  definitely  limit  the  amount  of  energy  the 
vapor  can  transport  down  the  pipe,  thus  limiting  the  opera- 
tion of  the  pipe.  The  heat  transfer  mechanism  employed  in 
the  rotating  heat  pipe  depends  on  the  latent  energy  of  the 
working  fluid.  The  heat  transfer  rate  in  the  pipe  may 
then  be  expressed  as: 

= "v  '’fg 
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where 


= mass  flow  rate  of  the  vapor,  Ibm/hr. 

From  continuity,  the  mass  flow  rate  of  the  vapor  may  be 
expressed  by: 


m 


= P. 


u 


V 


(2) 


where 


= velocity  of  the  vapor,  ft/sec 

2 

A = cross  section  area  for  the  flow,  ft  . 

Then 

q.  = p u A h,  . (3) 

^t  v v fg  ' ' 

The  vapor  velocity  is  taken  as  the  sonic  velocity,  so 
that : 

= c = /g^  kRT  (4) 

where 

c = sonic  velocity,  ft/sec 
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k 

R 

T 


acceleration  of  gravity,  32,1739  ■ 2 

Ibf-sec^ 

the  ratio  of  specific  heats, 
gas  constant,  f t-lbf/lbm-deg  R 
absolute  temperature,  deg  R 


Using  equations  (3)  and  (4) , the  sonic  limit  of  a 
typical  rotation  heat  pipe  can  be  computed,  as  shown  in 
Figure  2. 

2 , Boiling  Limit 

Kutateladze  [2]  has  postulated  that  the  transition 
from  nucleate  boiling  to  film  boiling  is  purely  a hydro- 
dynamic phenomenon.  Based  on  this  assumption  he  has 
obtained  the  correlation  formula  for  predicting  the  burnout 
heat  flux: 

^ 1/4 

‘^t  = ^b  ^fg  - Pv’  ^ 

where 

2 

= heat  transfer  area  in  the  boiler,  ft 
= density  of  the  vapor,  Ibm/ft^ 
h^g  = latent  heat  of  vaporization,  BTU/lbm 
o = surface  tension,  Ibf/ft 
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g 


2 

= acceleration  of  gravity,  ft/hr 
0^  = density  of  fluid,  Ibm/ft^ 

K = constant  value. 

Kutateladze ' s experimental  data  indicate  a value  of 
K between  0.13  and  0.19.  Using  equation  (5),  the  boiling 
limit  of  a typical  rotating  heat  pipe  can  be  computed,  as 
shown  in  Figure  2. 

3 . Entrainment  Limit 

Sakhuja  [3]  studied  the  flooding  constraint  in 
wickless  heat  pipes,  where  he  assumed  that: 

a)  the  vapor  and  liquid  are  in  counterflow  motion 
and  are  exposed  to  each  other; 

b)  the  falling  film  of  liquid  (under  gravity)  is 
subjected  to  a shear  stress  by  counterflowing 
vapor,  which  tends  to  retard  the  film; 

c)  the  liquid  film  remains  smooth  and  stable,  and 
the  shear  stress  is  usually  small; 

d)  the  buoyancy  forces,  due  to  density  differences 
between  the  vapor  and  the  liquid,  are  responsible 
for  maintaining  the  counterflow;  and 

e)  the  buoyancy  forces  are  balanced  with  dissipative 
effects  which  are  proportional  to  momentum 
fluxes  of  vapor  and  liquid  streams. 

He  has  obtained  a correlation  for  flooding: 


^x  hg^/gP(Pf -p^)p^ 

[1  + 


(6) 
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where 


D 

g 

P. 


V 


fg 


= heat  transfer  rate,  BTU/hr 

2 

= flow  area,  ft 

= dimensionless  constant,  0.725  for  tube 
with  sharp  edged  flange 

= inside  diameter  of  the  heat  pipe,  ft 

2 

= acceleration  due  to  gravity,  ft/hr 
= density  of  the  fluid,  Ibm/ft^ 

= density  of  the  vapor,  Ibm/ft^ 

= latent  heat  of  vaporization,  BTU/lbm. 


Using  equation  (6),  the  entrainment  limit  of  a 
typical  rotating  heat  pipe  can  be  computed  and  is  shown 
in  Figure  2. 

4 , Condensing  Limit 

The  solution  for  condensation  in  the  rotating  heat 
pipe  has  been  carried  out  by  Ballback.  [1].  He  modeled  the 
condenser  section  of  the  heat  pipe  as  a rotating,  truncated 
cone,  and  an  equation  to  compute  the  condensing  limit  was 
derived : 


% * 


. kp  Pp  h,  [T  -T 
2 f fg"  s w 1 


sin 


{ [R^  + L sin4i] 


8/3 


- R 


8/3  3/4 

} 


(7) 


where 

q^  = total  heat  transfer  rate,  BTU/hr 

k,  = thermal  conductivity  of  the  condensate 
film,  BTU/hr-f t-deg  F 
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0) 


L 


density  of  fluid,  Ibm/ft^ 
angular  velocity,  1/hr 

latent  heat  of  varporization , BTU/lbm 

saturation  temperature,  deg  F 

inside  wall  temperature,  deg  F 

viscosity  of  fluid,  Ibm/ft-hr 

half  cone  angle,  degrees 

minimum  wall  radius  in  the  condenser 
section,  ft 

length  along  the  wall  of  the  condenser,  ft. 


Equation  (7)  is  a function  of  the  physical  proper- 
ties of  the  working  fluid,  the  geometry  of  the  rotating 
heat  pipe  and  the  operating  speed.  Using  this  equation, 
the  condensing  limit  for  a typical  rotating  heat  pipe  can 
be  computed.  This  result  is  also  shown  in  Figure  2. 

It  is  clear  from  the  results  of  Figure  2,  that  for 
the  geometry  and  working  fluid  chosen,  and  under  normal 
circumstances,  the  heat  transfer  rate  is  limited  primarily 
by  the  amount  of  heat  which  can  be  transferred  through  the 
condenser  wall.  However,  other  limitations  may  become 
important  as  heat  pipe  geometry  and  operating  conditions 
are  altered. 


C.  THEORETICAL  INVESTIGATIONS  OF  FILM  CONDENSATION 
ON  SMOOTH  ROTATING  CONDENSER  SURFACES 

Sparrow  and  Gregg  [4]  were  the  first  to  study  condensa- 
tion on  a rotating  surface.  They  studied  condensation  on  a 
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f 


rotating  disk  at  uniform  surface  temperature,  T^,  situated 
in  a large  quiescent  body  of  pure  saturated  vapor.  They 
assumed  that  the  condensed  liquid  formed  a continuous  film 
on  the  disk  and  that  fluid  in  this  film  moved  radially 
outward  along  the  disk  under  the  action  of  the  centrifugal 
force  field.  They  arrived  at  the  following  solution  for 
the  situation  where  energy  convection  and  acceleration  terms 
are  negligible: 


where 


CO  ' 


= 0.904( 


c at/h;  ' 

p fg 
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= local  heat  transfer  coefficient, 

BTU/hr-f t^-deg  F 

= thermal  conductivity  of  condensate, 

BTU/hr-f t-deg  F 

2 

= kinematic  viscosity  of  condensate,  ft  /hr 

= angular  velocity,  1/hr 

= Prandtl  number,  -C  /k, 

P ^ 

= absolute  viscosity  of  condensate,  Ibm/ft-hr 
= specific  heat  of  condensate,  BTU/lbm-deg  F 
= disk  surface  temperature,  deg  F 
= saturation  temperature  of  vapor,  deg  F 
= T - T , deg  F 

= latent  heat  of  vaporization,  BTU/lbm. 
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Equation  (8)  predicts  that  the  local  heat  transfer 
coefficient  is  constant  across  the  disk  and  does  not  vary 
with  radial  position. 

Sparrow  and  Hartnett  [5]  continued  the  research 
using  a rotating  cone  with  the  same  assumptions  used  earlier. 
Their  result  for  the  heat  transfer  coefficient  was: 


h /hj . , 

cone  disk 


1/2 


(sin  't>) 


(9) 


where  o is  the  internal  half 'angle  of  the  cone.  As  pointed 
out  by  the  authors,  however,  equation  (9)  applies  only  to 
cones  that  are  not  too  slender.  Using  their  result,  the 
mean  Nusselt  number  for  a cone  can  be  written  as: 


Nu 
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fg 
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Leppert  and  Nimmo  [6,7]  studied  laminar  film 
condensation  on  a surface  normal  to  body  or  inertial  forces 
such  as  a finite  horizontal  surface,  or  a rotating  drum. 

They  found  that  the  local  heat  transfer  coefficient  varies 
with  position,  and  the  mean  Nusselt  number,  for  a constant 
surface  temperature,  with  negligible  surface  tension  effects, 
is : 


Nu  = 0.64  Sh^'^^ 

m 


(11) 
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where 


Nu  = mean  Nusselt  number. 


h L 
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2 V,  r 3 

g P L 

Sh  = Sherwood  number,  — , — 

k AT  u 


Equation  (11)  can  be  applied  to  the  case  of 

condensation  on  the  inside  of  a rotating  cylinder  of  radius, 

2 

r,  where  g is  replaced  by  rw  . 

The  first  theoretical  investigation  conducted  at 
the  Naval  Postgraduate  School  was  performed  by  Ballback 
[1].  His  analysis  was  a Nusselt  type  analysis  for  film 
condensation  on  the  inside  surface  of  a rotating  truncated 
cone  of  half  angle  cone,  <p , with  a geometry  defined  in 
Figure  3.  His  assumptions  were  similar  to  those  used  in 
the  Nusselt  analysis  of  film  condensation  on  a vertical 
wall,  namely: 

1)  steady  state  condition, 

2)  film  condensation,  not  dropwise  condensation, 

3)  laminar  flow  of  condensate, 

4)  constant  fluid  properties, 

5)  no  condensate  subcooling, 

6)  no  momentum  change  through  the  condensate, 

7)  no  liquid-vapor  interfacial  shear  stresses, 

8)  linear  distribution  of  temperature  through  the 
condensate  film, 

9)  negligible  force  of  gravity,  and 

10)  constant  wall  temperature. 
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Additionally  Ballback  assumed: 


11) 

a very  small  film  thickness,  5 (x) , 

12) 

no  radial  distribution  of  pressure 
vapor  space. 

in  the 

13) 

negligible  density  of  the  vapor. 

14) 

negligible  velocity  gradients  in  the 
circumferential  direction,  and 

15) 

very  small  slope  of  the  condensate 
d6/dx  when  compared  with  tan  1>. 

film. 

He  used  the  boundary  condition,  that  5 (x)  must  be 
zero  at  X = 0,  and  arrived  at  the  following  solution  for 
the  local  film  thickness: 


5 (x)  = 


Uf(T^-T„) 


•[1  - { 


sin  ;()  h 


fg 


R + X sin  (t 
o ^ 


8/3  1/4 

} ]} 


(12) 


where 


>5  (x) 


<P 

R 

o 


film  thickness,  ft 

thermal  conductivity  of  the  fluid, 
BTU/hr-f t-deg  F 

viscosity  of  the  fluid,  Ibm/ft-hr 

density  of  fluid,  Ibm/ft^ 

angular  velocity,  1/hr 

inside  wall  temperature,  deg  F 

saturation  temperature,  deg  F 

latent  heat  of  vaporization,  BTU/lbm 

half  cone  angle,  degrees 

minimum  wall  radius  in  the  condenser 
section,  ft. 
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The  local  heat  transfer  coefficient,  h(x),  is 


given  by; 


h (x)  = 


(x) 


3 2 2 

^ p (jj  sin  <t>  h^ 

{£  ^ f . fg 

3 U,(T^-TJ 


1/4 


- 1 - [ 


+ X sin  t> 


8/3  ^ 

] } 

(13) 


The  mean  heat  transfer  coefficient  over  the  whole 
surface  area  of  the  truncated  cone  is  given  by: 


m 


ilTTR  +Lsin  pi  ''  2-rh{x)[R  +x  sin  p]  dx 

o 0 


(14) 


where 

L = length  along  the  wall  of  the  condenser,  ft. 

Integrating  equation  (14)  yields  a value  for  the 
mean  Nusselt  number: 
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o _ L sin  () 

® ~ ~R ' 

o 

Equation  (15)  can  be  put  into  the  form: 

Nu^  = (Nu^Uh  (17) 

m m S H 

where  (Nu  ) _„  is  the  result  of  Sparrow  and  Hartnett  [5] 
in  on 

for  a cone,  as  given  by  equation  (10).  The  above  result 
was  confirmed  as  a special  case  of  the  later  work  of  Dhir 
and  Lienhard  [8],  who  studied  laminar  film  condensation  on 
axisymmetric  bodies  in  nonuniform  gravity  fields. 

Daley  [9]  was  able  to  improve  on  the  theory  intro- 
duced by  Ballback  by  consideration  of  the  thermal  resis- 
tances in  the  condenser  wall  and  on  the  outside  surface. 
Both  Ballback  and  Daley,  however,  neglected  the  drag  due 
to  vapor  friction  on  the  counterflowing  condensate  film. 

Newton  [10]  extended  the  analysis  performed  by 
Ballback  and  Daley  to  include  the  effects  of  vapor  pressure 
drop  and  interfacial  shear  stress  between  vapor  and  conden- 
sate. Using  the  same  geometry  in  Figure  3,  he  assumed  that 
the  momentum  changes  within  the  condensate  film  were  small 
and  arrived  at  the  following  differential  equations  for 
condensate  film  thickness,  vapor  velocity  and  heat  transfer 
rate : 
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where 


P^(x) 


P2  (x) 


v(x) 

T_ 


t 

5 (x) 
k^ 


6^x)  _ 4 

(R^  + xsin  4)  ) — - -jjcos  4>  S (x)  , 

^ 3 

(R  + xsin  <p  )-^4^  “ cos  <p  6 (x)  , 

° 3 

velocity  of  vapor,  ft/hr 

saturation  temperature  of  vapor,  deg  F 

ambient  coolant  temperature,  deg  F 

thickness  of  the  condenser  wall,  ft 

film  thickness  of  the  condensate,  ft 

thermal  conductivity  of  condensate  film, 
BTU/hr-ft-deg  F 
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R(x) 
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thermal  conductivity  of  condenser  wall, 
BTU/hr-f t-deg  F 

external  heat  transfer  coefficient, 
BTU/hr-f t^-deg  F 

latent  heat  of  vaporization,  BTU/lbm 
friction  factor,  dimensionless, 
vapor  Reynolds  number,  2R (x) p^v (x) 

radius  along  the  condenser  wall, 

R + xsin  p - 5cos  (p 
o 

half  cone  angle,  degrees 
density  of  vapor,  Ibm/ft^ 
viscosity  of  vapor,  Ibm/ft-hr. 


The  integrations  were  started  by  assuming  initial 

conditions  5=5.  and  v = 0 at  x = 0.  With  these  initial 

1 

conditions  ^ — = tan  i at  x = 0.  Newton  obtained  results 
dx 

for  half  cone  angles  of  0,  0.1,  0.2  and  0.3  degrees.  As 
the  half  cone  angle  increased,  the  solution  of  the  differ- 
ential equations  became  extremely  sensitive  to  the  initial 
values  of  5^.  He  also  noted  that  his  results  for  5 (x) 

showed  that  (except  near  the  ends  of  the  condenser) 
d 6 

^ <<  tan  p.  He,  therefore,  made  this  approxim.ation  and 
simplified  the  model  further  using  an  order  of  magnitude 
analysis  of  the  terms  in  equations  (18)  and  (19).  With 
these  simplifications,  equations  (18)  and  (19)  were 
reduced  to: 
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4 


2 

PpR(x) 

- o^v(x)-2-^ = 0 (22) 


Equation  (21)  was  numerically  integrated  using  the 
initial  condition  that  5^  = 0 and  v = 0 at  x = 0 and  using 
the  cubic  equation  for  6(x),  equation  (22).  After  the 
integration  was  completed  for  5 (x)  and  v(x),  the  energy 
equation,  equation  (20)  , was  integrated  over  the  length  of 
the  condenser  to  find  the  total  heat  transfer  rate  of  the 
heat  pipe.  This  theory  was  compared  to  the  experimental 
results  of  Schafer  [15]  and  Tucker  [20],  predicting  perfor- 
mance to  within  i20  percent  for  water,  ethyl  alcohol  and 
Freon  113. 

Daniels  and  Al-Jumaily  [11,12]  extended  the  analysis 
of  Ballback  by  including  drag  effects  due  to  the  vapor 
friction.  Their  solution  for  t.he  local  film  thickness 
distribution  along  the  condenser  is  given  by: 
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Sherwood  number, 
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drag  friction  number. 


OT  h^  X cos  t 
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two  phase  Reynolds  number,  — ~ — 

liquid  density,  Ibm/ft^ 

film  thickness,  ft 

angular  velocity,  1/hr 

half  cone  angle,  degrees 

shear  stress  at  liquid  vapor  interface, 
Ibf/ft^ 


u 


liquid  condensate  velocity,  ft/hr 

2 

acceleration  of  gravity,  ft/hr 
latent  heat  of  vaporization,  BTU/lbm. 


At  any  section  x in  the  condenser,  since  heat  is 
transferred  through  the  condensate  layer  by  conduction 
only,  the  local  heat  flux  = kAT/6 (x)  and  equation  (23) 
can  be  re-arranged  as  follows: 


q4  , iRe  (“I)  2q2  ^ Ip,  3q  . sh  (^AI)  ^ » 0 (24 

X 4 VX  X X 3 X X X X X 
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where 


= local  heat  flux,  BTU/hr-ft^. 

This  equation  relates  the  three  main  parameters 

controlling  the  rotating  heat  pipe  operation  assuming 

fixed  geometry  and  working  fluid.  The  mean  Nusselt  number, 

Nu  , for  this  analysis  was, 
m 


where  suffix  L refers  to  quantities  evaluated  at  the  end 
of  the  condenser  section  x = L. 

Marto  [13]  compared  the  heat  transfer  rates  of 
Daniels  and  Al-Jumaily  with  Ballback  and  the  approximate 
numerical  solution  to  equations  (20)  , (21)  and  (22)  . For 

the  case  of  Freon  113,  equations  (20),  (21)  and  (22), 

which  include  both  effects  of  interfacial  shear  as  well  as 
vapor  pressure  drop,  gave  slightly  lower  heat  transfer 
rates  than  the  analysis  of  Daniels  and  Al-Juraaily  which 
includes  only  interfacial  effects.  Also  the  analysis  of 
Ballback  predicts  slightly  larger  results  than  the  computer 
analysis  but  lower  results  than  those  of  Daniels  and 
Al-Jumaily.  This  may  be  due  to  differences  in  the  Freon  113 
fluid  properties  used  or  due  to  the  geometrical  inconsistency 
in  matching  all  the  geometrical  dimensions  as  stated  in  the 
analysis  of  Daniels  and  Al-Jumaily.  For  the  case  of  water, 


Ballback's  results  again  predict  a slightly  higher  heat 
transfer  rate  than  the  approximate  numerical  solution,  due 
to  his  neglecting  effects  of  interfacial  shear  and  vapor 
pressure  drop.  When  these  effects  are  neglected  in  the 
approximate  numerical  solution,  the  results  agree  very 
well  with  those  of  Ballback.  For  these  reasons,  the 
effects  of  interfacial  shear  and  vapor  pressure  drop  can 
safely  be  neglected. 

Roetzel  [14]  studied  condensation  in  steam  heated 
fast  rotating  paper  drying  drums,  as  shown  in  Figure  4. 

In  his  analysis,  the  heat  transfer  is  determined  mainly  by 
the  method  of  condensate  removal.  The  condensate  is  collected 
in  a circumferential  groove  on  the  inside  surface,  and  from 
there  it  is  removed  by  means  of  a stationary  syphon  with 
its  tip  in  the  groove.  The  mean  film  thickness  is  slightly 
smaller  than  with  a rotating  syphon  because  the  condensate 
is  sucked  off  from  a lower  level.  The  solution  for  the 
case  of  a drum  with  an  infinite  radius  of  curvature  (The 
inside  surface  is  exactly  cylindrical.)  with  a uniform 
temperature  difference  is  given  by; 

1/5 

Nu^  = (0.86  E)  (26) 


where 


Nu  = Nusselt  number  for  the  cylindrical  drum. 
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= thermal  conductivity  of  condensate, 
BTU/hr-ft-deg  F 

= Sherwood  number 

y 

= acceleration  normal  to  the  wall,  ft/hr 
= density  of  condensate,  Ibm/ft^ 

= latent  heat  of  vaporization,  BTU/lbm 
= total  length  of  condensate,  ft 
= kinematic  viscosity  of  condensate,  ft  /hr, 


This  equation  can  be  written  in  terms  of  the 
Sherwood  number: 


Nu  = 0.678  Sh 


1/5 


(27) 


and  compares  favorably  to  the  analysis  of  Nimmo  and  Leppert, 
equation  (11) . 

For  the  case  of  finite  radius  of  curvature,  Roetzel's 
solution  becomes: 


= 

m 


23/16 
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(0.86E)23/20j4/23 


(28) 
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where 


R = radius  of  curvature,  ft,  and 
L = total  flow  length  of  condenser,  ft. 


This  equation  can  be  put  in  terms  of  Sh,  Nu^,  R,  and  L as: 


Nu^  = NuJ1.  1.924(|)23/16  3^23/80^4/23  ^ (29) 


Thus  the  Nusselt  number  can  be  directly  determined  when 
the  ratio  L/R  and  the  value  of  the  dimensionless  group,  E, 
are  known.  For  the  truncated  cone,  rotating  heat  pipe,  as 
shown  in  Figure  3,  equation  (29)  can  be  applied  to  compute 
the  heat  transfer  by  defining  the  radius  of  curvature 
R = L/tan  0 ; the  solution  then  becomes: 

Nu^  = Nu_^  [1  + 1.924  tan  $ (30) 


where 


h = half  cone  angle,  degrees 

Nu^  = Nusselt  number  of  the  cylindrical  drum, 
equation  (11) . 

Equation  (30)  has  never  been  compared  with  experi- 
mental data,  so  using  the  existing  rotating  heat  pipe  data 
together  with  new  experimental  data  at  different  half  cone 
angles  will  provide  a check  on  its  usefulness. 
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These  theoretical  investigations  are  summarized 
below : 

1.  The  vapor  pressure  drop  and  interfacial  shear 
stress  effects  between  vapor  and  condensate  can 
be  neglected. 

2.  The  analysis  of  Roetzel  [14]  appears  to  be  the 
most  appropriate  for  use  in  a rotating  heat  pipe. 

3.  The  important  parameters  are: 

(a)  speed, 

(b)  geometry,  and 

(c)  Sherwood  number,  which  includes  the  properties 
of  the  working  fluid. 

D.  ANALYSIS  OF  FINNED  ROTATING  CONDENSER  SURFACES 

The  use  of  axial  fins  on  the  inner  surface  of  the  con- 
denser can  improve  heat  transfer  performance  through  two 
mechanisms:  a decrease  in  condensate  film  thickness  near 

the  apex  of  each  fin  and  an  increase  in  inner  wall  surface 
area.  Schafer  [15]  developed  an  analytical  model  of  the 
internally  finned  condenser  for  the  case  of  a fin  of  tri- 
angular profile  as  shown  in  Figure  5,  with  a specific 
geometry  as  shown  in  Figure  6.  In  the  analysis,  he  assumed 
one-dimensional  heat  conduction  from  the  inside  to  the 
outside  surface  and  a linearly  decreasing  surface  tempera- 
ture profile  from  the  fin  apex  to  the  midpoint  of  the 
trough.  Results  from  his  analysis  indicated  significantly 
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improved  heat  transfer  rates  for  the  finned  condenser 
compared  to  the  smooth  condenser. 

Corley  [16]  developed  a numerical  solution  for  this 
case  using  the  Finite  Element  method  to  account  for  two- 
dimensional  conduction  effects  within  the  wall  for  the 
triangular  fin  profile  and  a parabolic  wall  surface  tempera- 
ture distribution  along  the  fins.  His  results  indicated 
significantly  improved  heat  transfer  performance  with  heat 
transfer  rates  about  75  percent  above  those  predicted  by 
the  one-dimensional  linear  temperature  distribution  analy- 
sis of  Schafer  [15].  However,  Corley  noted  that  an  error 
of  50  percent  in  the  calculation  of  apex  heat  transfer  rate 
might  then  give  a total  heat  transfer  error  of  as  high  as 
15  percent.  Two  possible  solutions  to  the  boundary  con- 
dition problem  were  recommended.  The  first  was  to  increase 
the  number  of  Finite  Elements  in  the  analysis,  and  the 
second  solution  was  to  consider  a wall  geometry  as  shown 
in  Figure  7. 

E.  THESIS  OBJECTIVES 

The  objectives  of  this  thesis  were  therefore: 

1)  to  find  out  which  existing  theory  best  correlates 
Rotating  Heat  Pipe  data  by  comparing  these  theories 
with  existing  data  and  new  experimental  data  of 
this  thesis,  and 

2)  to  further  study  finned  condensers  by  improving 
the  Finite  Element  results  obtained  earlier. 


41 


err 


FINITE  ELEMENT  ANALYSIS  AND  RESULTS 


II . 

A.  REVIEW  OF  EARLIER  ANALYSES 

Schafer  [15]  studied  the  one-dimensional  heat  transfer 
solution,  and  Corley  [16]  studied  the  two-dimensional  heat 
transfer  solution  to  finned  condensers.  Both  used  the 
same  assumptions  and  boundary  conditions,  based  upon  the 
work  of  Ballback  [1],  which  are  similar  to  those  used  in 
the  Nusselt  analysis  of  film  condensation  on  a vertical 
wall.  The  more  important  of  these  assumptions  are  listed 
below ; 

1)  steady  state  operation, 

2)  film  condensation,  as  opposed  to  dropwise  condensation, 

3)  laminar  condensate  flow  along  both  the  fin  and 
the  trough, 

4)  static  balance  of  forces  within  the  condensate, 

5)  one-dimensional  conduction  heat  transfer  through 
the  film  thickness  (no  convective  heat  transer) , 

6)  no  liquid-vapor  interfacial  shear  stresses, 

7)  no  condensate  subcooling, 

8)  zero  heat  flux  boundary  conditions  on  both  sides 
of  the  wall  section,  as  shown  in  Figure  6, 

9)  saturation  temperature  at  the  fin  apex, 

10)  zero  film  thickness  at  the  fin  apex, 

11)  negligible  curvature  of  the  wall. 

The  complete  development  of  Schafer's  finned  condenser 
heat  transfer  theory  is  presented  in  his  thesis  [15].  His 
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results  are  used  for  comparison  purposes  in  this 
thesis . 

Corley's  study  [16]  considered  an  internal  fin  geometry 
as  shown  in  Figure  6.  Two-dimensional  conduction  through 
the  condenser  wall  was  considered  in  the  heat  transfer 
analysis  by  using  a Finite  Element  solution.  Instead  of 
using  a linearly  decreasing  temperature  profile  along  the 
wall  from  the  apex  of  the  fin  to  the  midpoint  of  the 
channel  as  Schafer  did,  a parabolic  distribution  was  assumed 
to  calculate  average  nodal  convective  coefficients.  The 
eight  node  finite  element  mesh  used  in  the  two-dimensional 
analysis  of  the  axially  finned  condenser  is  shown  in  Figure 
8.  A computer  program  for  the  case  of  two-dimensional  con- 
duction was  developed  using  a modified  version  of  the  three 
dimensional,  isoparametric  Finite  Element  program  assembled 
by  Lew  [ 17 ] . 

A parabolic  temperature  distribution  is  assumed  to  exist 
along  the  fin  surface  using  the  boundary  condition,  T = T^ 
at  z = 0,  since  5(0)  = 0.  Therefore; 

2 

T'z)  = a.z  +b,z+T 
w ' 1 Is 


or 


^s  - 


T (z)  = 

w 


-aiz 


- b^z 


(31) 
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+ X sin  J) 


t = 0.03125  in. 
b = 0.025  in. 

R = 0.762  in. 

0 = 1® 


Figure  8.  Condenser  Geometry  Considered  in 

Two-Dimensional  Analysis  with  2 Finite 
Elements 


where 


and  are  unknown  constants,  and 

z is  the  distance  along  fin  surface  from  apex  of 
fin,  ft. 

From  Schafer's  analysis  [15]  a differential  equation 
for  fin  condensate  thickness,  5 (z) , can  be  derived: 


6 (z)d6  (z) 


k.(T^,  - T (z)  ) dz 

I b W t 


2 2 

D-oj  (R  + xsin  l))h^  cos  3 cos  a 
f o fg 


(32) 


where 


6 (z) 
”^f 

^f 

^f 


00 


R 


O 


X 


z 


a 


the  fin  condensate  thickness,  ft 

thermal  conductivity  of  working  fluid, 
BTU/hr-ft-deg  F 

density  of  working  fluid,  Ibm/ft^ 

viscosity  of  working  fluid,  Ibm/ft-hr 

angular  velocity,  1/hr 

minimum  radius,  ft 

distance  along  apex  of  fin,  ft 

distance  along  the  fin,  from  the  apex,  ft. 

half  cone  angle,  degrees 

fin  half-angle,  degrees 

latent  heat  of  vaporization,  BTU/lbm. 
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Equation  (31)  is  then  substituted  in  (32),  and  the 
result  is  integrated  from  z = 0 to  z to  give  an  equation 
for  6 (z  ) : 


<5  (z) 


4kf(-a^  X - h V 


2 2 

p -u)  (R  + xsin  <t> ) h^ 
f o fg 


U.  1/4 

] 

cos  Ip  cos  a 


(33) 


The  local  convective  heat  transfer  coefficient  assuming 
one-dimensional  heat  conduction  through  the  thin  condensate 
film  is  then: 


h (zl 


f in 


’'f 
5 (z) 


= [ 


3 2 2 

(R  +xsinP)h^ 
f f o fg 


cos  p cos  a 1/4 
] 


4p j (-a^ 


- b 


1 T 


(34; 


which  is  applicable  from  z = 0 to  z*  = z^  - 5*/cos  a , 
where  5*  is  the  trough  condensate  thickness,  as  shown  in 
Figure  8.  Across  the  width  of  trough,  from  z = z^  to  the 
mid  channel  position,  the  heat  transfer  coefficient  is: 


trough  5 * ( x ) 


(35) 


where  6*(x)  is  considered  constant  across  the  trough  width. 
Based  upon  a check  of  the  wall  temperatures  generated  from 
Schafer's  computer  program,  the  fluid  temperature  difference 
across  the  width  of  the  trough  was  predicted  to  be  small. 
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Fluid  thermal  conductivity  values  were  assumed  to  remain 


relatively  constant.  As  a result,  h(x)^  , is  assumed 

trough 

to  be  constant  across  the  trough  width  and  is  calculated 
using  equation  (35)  and  a representative  average  inner 
wall  temperature  to  determine  The  convective  coeffi- 
cient along  the  outside  surface  of  the  wall,  h is  also 

out 

assumed  to  be  constant.  The  modified  trough  mass  flow 
equation  becomes : 


2 2 

0 w 


(R  + xsin  t>)  0 * (x)  sin  i) 


3 U. 


2 

[5*(x)t  +5*(x)tan  i ] 


(36) 


Within  the  computer  program,  the  heat  transfer  rate 
for  the  entire  condenser  is  determined  by  dividing  its 
axial  length  into  one  hundred  equal  length  increments  of 
length.  Ax  (See  Figure  6.),  solving  for  the  incremental 
heat  transfer  rate  using  an  iterative  procedure,  and 
summing  these  incremental  rates  to  yield  a total  heat 
transfer  rate.  A flow  chart  description  of  the  computer 
program  used  in  Corley's  analysis  is  given  in  Figure  9. 

At  each  increment,  constant  values  of  a^^  and  b^^  are  deter- 
mined using  a Lagrangian  interpolation  fit  of  the  tempera- 
tures, at  nodes  1,  3 and  6 as  calculated  from  the  previous 
iteration.  For  the  first  increment,  a^  and  b^  are  calcu- 
lated from  an  initial  guess  of  the  temperatures  at  these 
nodes.  Convective  coefficients  at  several  locations  along 
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* 


Figure  9. 


computer  Program  Flowchart  of  ^o-Dimensional 
Conduction  Analysis  of  Corley  lloj 
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the  fin  are  determined  using  equation  (34).  Average 


nodal  values,  ^3avg'  average  h of  the 

left-hand  boundary  region  associated  with  node  6,  in  Figure 
8,  are  calculated  using  the  position  dependent  values  of 
h(z)^^^  in  a Simpson's  Rule  approximation.  In  the  trough, 
is  calculated  from  equation  (35)  using  the  5* 
value  calculated  from  the  previous  increment.  This  value 
is  used  as  the  average  convective  coefficient  for  nodes  9 
and  12  and  the  right  hand  boundary  region  of  node  6.  The 
average  convective  coefficients  calculated  for  the  left- 
and  right-hand  boundary  regions  of  node  6 are  weighted 
by  their  respective  region  lengths.  The  two  values  are 
summed  to  yield  the  average  convective  coefficient  for  node 
6,  remaining  boundary  conditions  are  applied, 

and  the  incremental  heat  transfer  rate  per  unit  condenser 
length,  , and  nodal  temperature  values  are  generated 
from  the  finite  element  subroutines.  If  the  resultant 
heat  transfer  rate  converges  to  within  the  criterion  of 
0.05  percent  change  from  t.he  value  of  the  previous  iteration. 


iQ-i  - Qi-i' 

-3 <<  0.0005  , 

I ^3  ! 


(37) 


then  the  incremental  problem  is  deemed  solved  and  is 
set  equal  to  the  final  iteration,  Q ^ . If,  on  the  other 
hand,  the  convergence  criteria  is  not  met,  newly  calculated 
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nodal  temperature  values  are  used  to  establish  a new  set 
of  inner  wall  convective  coefficients,  and  the  finite  element 
solution  is  repeated  until  the  heat  transfer  rate  converges 
to  within  the  acceptable  0.05  percent  tolerance.  Because 
Qj  is  initialized  at  zero  for  each  increment,  the  iteration 
procedure  is  repeated  at  least  twice  for  every  increment. 
Incremental  mass  flow  rate  is  calculated  using: 


AM 


tot 


2 ^ Ax 

hf: 

fg 


(38) 


This  value  is  added  to  the  mass  flow  rate  from  previous 
increments.  The  subsequent  interval's  trough  condensate 
thickness,  5*(x),  is  calculated  using  a polynomial  root- 
finder  subroutine  to  solve  equation  (36) . This  incremen- 
tation is  continued  until  the  entire  length  of  the  conden- 
ser is  traversed.  Incremental  heat  transfer  rates  are 
summed,  and  a total  heat  transfer  rate  becomes: 


'tot 


100 

= 2 N_.  E (Q-  Ax) 


(39) 


where  N, • is  the  number  of  fins  along  the  inner  wall 
circumference.  To  start  the  incrementation,  two  initial 
assumptions  are  made.  The  initial  value  of  5*  = 5|  is 
taken  from  the  analysis  by  Sparrow  and  Gregg  of  condensa- 
tion on  a rotating  disk  [4] . In  addition,  at  the  first 
increment,  initial  guesses  of  fin  temperatures  at  nodes  1, 
3 and  6 are  made  in  order  to  begin  the  solution. 
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B.  MODIFICATION  TO  THE  FINITE  ELEMENT  ANALYSIS 

This  thesis  is  a continuation  of  Corley's  [16]  study 
of  the  addition  of  internal  axial  fins  to  the  condenser 
surface  of  a rotating  heat  pipe.  He  recommended  that  one 
way  to  approximate  more  closely  the  actual  heat  transfer 
rates  would  be  by  increasing  the  number  of  elements  for 
the  same  geometry.  Since  Corley  used  one  hundred  equal 
increments  of  length  to  compute  the  heat  transfer  rates, 
the  first  step  in  this  thesis  was  to  reduce  the  number  of 
increments  to  fifty  and  twenty  five,  to  minimize  the  com- 
puter time.  Using  Corley's  program,  at  a speed  of  3000 
RPM,  with  external  convective  coefficients  of  1000,  1500 
and  5000  BTU/hr-f t^-deg  F and  fin  half  angles  of  10,  30 
and  50  degrees,  the  results  for  water  as  the  working 
fluid  are  shown  in  Tables  2,  3 and  4.  It  is  clear  that 
reducing  the  number  of  increments  to  twenty  five  had  little 
effect  upon  the  heat  transfer  results  (approximately  3 
percent  error) . 

To  reduce  any  error  in  the  calculation  of  the  heat 
transfer  rate  at  the  apex  of  the  fin,  the  number  of  Finite 
Elements  has  to  be  increased  between  the  apex  of  the  fin 
and  the  trough,  as  shown  in  Figure  10.  The  geometry  becomes 
3 elements  and  18  nodes.  The  computer  program  was  therefore 
modified  to  solve  this  case,  using  the  same  boundary  condi- 
tions and  assumptions.  These  results  are  shown  in  Tables  5 
through  10.  Upon  dividing  element  in  Figure  10  into 
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TABLE  2 


Theoretical  Finned  Condenser  Heat  Transfer 
Rate  Versus  Number  of  Axial  Increments 


h ^ = 1000 
out 

BTU/hr-ft^-deg 

F 

At  3000 

RPM 

'in  Half 

No.  of  T ^ 

_ . sat 

Fins 

No.  < 

af  Axial 

Increments 

Angle 

25 

50 

100 

10 

276  100 

10273.1 

10039.0 

9922.5 

150 

27355.8 

26732.6 

26422.4 

200 

44427.2 

43415.0 

42911.2 

30 

84  100 

10095.3 

9865.3 

9750.3 

150 

26866.0 

26253.8 

25949.1 

200 

43622.0 

42627.9 

42133.2 

50 

40  100 

9870.9 

9643.5 

9530.3 

150 

26255.7 

25650.8 

25349.9 

200 


42623,3 


41641.2  41152.6 


TABLE  3 


Heat  Transfer 
Increments 


Theoretical  Finned  Condenser 
Rate  Versus  Number  of  Axial 


h ^ = 1500  BTU/hr 
out 

-f t^-deg 

F 

At  3000 

RPM 

Fin  Half  No.  of 

Angle  Fins 

sat 

No. 

25 

of  Axial 
50 

Increments 

100 

10  276 

100 

15195.2 

14849.2 

14677.1 

150 

40427.4 

39507.2 

39049.2 

200 

65632.8 

64138.9 

63395.3 

30  84 

100 

14808.4 

14471.6 

14304.0 

150 

39368.6 

38473.3 

38027.8 

200 

63896.8 

62443.8 

61720.6 

50  40 

100 

14368.5 

14038.7 

13874.6 

150 

38174.4 

37298. 1 

36862.  1 

200 

61943.9 

60521.6 

59814.1 
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TABLE  4 

Theoretical  Finned  Condenser  Heat  Transfer 
Rate  Versus  Number  of  Axial  Increments 

= 5000  BTU/hr-ft^-deg  F At  3000  RPM 


’in  Half 

No . of 

No 

1.  of  Axial 

Increments 

Angle 

Fins 

25 

50 

100 

10 

276 

100 

45685.4 

44652.3 

44138.1 

150 

120807.7 

118075.1 

116715.1 

200 

195703.4 

191270.2 

189066.4 

30 

84 

100 

42327.3 

41377.0 

40904.1 

150 

111871.5 

109356.8 

108106.0 

200 

181198.5 

177121.2 

175092.6 

50 

40 

100 

39287.1 

38402.4 

37962.4 

150 

103724.0 

101386.0 

100222.1 

200 

167950.6 

164155.9 

162268.8 
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Comparison  of  Theoretical  Finned  Condenser 
Heat  Transfer  Rate  at  300()  RPM,  for 
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Comparison  of  Theoretical  Finned  Condenser 
Heat  Transfer  Rate  at  1000  RPM,  for 
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Comparison  of  Theoretical  Finned  Condenser 
Heat  Transfer  Rate  at  3000  RPM,  for 
h = 5000  BTU/hr-f t2-deg  F 


two  equal  parts,  as  shown  in  Figure  11,  the  geometry 
becomes  4 elements  and  23  nodes.  The  computer  program 
was  modified  again  to  solve  this  case,  using  the  same 
boundary  conditions  and  assumptions  as  before.  Results 
for  this  case  are  also  tabulated  in  Tables  5 through  10. 

C.  DISCUSSION  OF  NUMERICAL  RESULTS 

Computer  results  for  the  finned  copper  condenser  were 
generated  for  three  outside  convective  coefficients,  three 
fin  angles  and  two  rotation  rates.  The  working  fluid  used 
throughout  the  analysis  was  water.  Tables  5 through  10 
compare  the  results  of  using  2,  3 and  4 elements,  with  those 
of  Schafer  [15 ] . 

Using  the  Finite  Element  solution,  the  heat  transfer 
rates  change  less  than  10  percent,  as  the  number  of  elements 
change,  so  the  heat  transfer  rates  in  the  finned  condenser 
can  be  approximated  by  2 elements-,  and  the  error  in  the 
calculation  of  the  heat  transfer  rates  at  the  fin  apex  can 
be  neglected.  In  comparing  the  smooth  condenser  to  the 
finned  condenser,  the  most  dramatic  improvement  occurs 
for  a = 50  degrees  with  40  fins.  It  is  clear  in  this 
situation  that  Element'^'^smooth  greater  than 

^4  Element^^smooth ■ 

The  heat  transfer  rate  was  found  to  vary  slightly  with 
rotational  speed  in  a finned  condenser.  This  result  agrees 
with  the  results  obtained  earlier  by  Schafer  [15].  Schafer 
suggested  that  the  invariance  in  heat  transfer  rate  with 
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the  rotation  rate  was  due  to  the  cancelling  effects  of 

decreasing  film  thickness  and  decreasing  temperature 

difference  across  the  film  thickness  as  rotation  rate  is 

increased.  The  influence  of  external  heat  transfer 

coefficient,  is  shown  in  Tables  5 through  10.  A 

large  external  heat  transfer  coefficient  (h  ^ = 5000 
^ out 

2 

BTU/hr-ft  -deg  F)  is  seen  to  produce  the  largest  heat 
transfer  rates. 
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Ill . 


EXPERIMENTAL  PROGRAM 


A.  EXPERIMENTAL  EQUIPMENT 

1 . Description  of  Equipment 

The  equipment  used  to  support  the  experimental 
program  was  the  same  equipment  used  by  Loynes  [18]  and 
Wagenseil  [19].  Figures  12  and  13  depict  the  experimental 
apparatus.  The  evaporator  was  a copper  cylinder  with 
0-ring  seals  at  each  end.  Heat  input  to  the  evaporator 
was  provided  by  a resistance  heating  element  wrapped 
around  the  evaporator.  The  power  supply  for  the  heater 
element  was  changed  from  load  bank  resistors  to  a solid 
state  phase  amplifier  power  controller  to  reduce  the  space 
and  inconvenience  of  the  load  banks.  The  model  PA-1-2490 
phase/amp  is  a single  phase  silicon  controlled  rectifier 
power  amplifier  designed  to  control  and  regulate  A-C  power. 

A six-volt  battery  was  used  as  a power  supply  for  the 
controller  circuit. 

The  condenser  was  surrounded  by  a sealed  insulated 
box.  Gar lock  seals  at  each  end  of  the  condenser  kept  the 
cooling  water  from  leaking  out.  The  cooling  water  system, 
as  shown  schematically  in  Figure  14,  was  four  supply  tubes, 
0.25  inch  diameter  with  0.031  inch  holes  drilled  every 
0.20  inch.  The  supply  tubes  ran  the  length  of  the  conden- 
ser, parallel  to  the  axis.  Filtered  and  softened  tap  water 
was  used  as  the  cooling  fluid.  The  inlet  and  outlet  cooling 
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Figure  13.  Schematic  Diagram  of  Heat  Pipe  Apparatus 


Cooling  water 
inlet 


Figure  14.  Impingement  Cooling  Schematic 


water  temperatures  were  measured  along  with  the  cooling 
water  flow  rate.  The  heat  pipe  was  rotated  in  two  bearings 
by  a three  phase  variable  speed  motor  with  a V-belt  drive. 

A small  drop-oiler  was  installed  over  the  main 
bearing  housing  to  lubricate  the  main  bearing.  This  unit 
consisted  of  a small  reservoir  with  a needle  valve  for  flow 
regulation.  A sight  glass  permitted  observation  of  the  oil 
flow  from  the  reservoir  to  the  main  bearing.  A flow  rate 
of  one  drop  every  20-30  seconds  was  found  to  be  adequate. 

The  end  of  the  evaporator  opposite  the  condenser  contained 
a viewport  consisting  of  two  glass  disks  with  an  air  space 
between  them  to  prevent  fogging.  A strobe  light  shined 
through  this  viewport  and,  when  synchronized  to  the  rota- 
tional speed,  allowed  observation  of  the  boiling  action  in 
the  fluid  annulus  and  the  condensate  in  the  condenser. 

2 . I.nstrumentation 

The  condenser  wall  was  grooved  out  externally  in 
nine  places  in  Figures  15  and  16.  The  grooves  were  spaced 
40  degrees  apart  and  had  a depth  of  0.010  inch.  Welded 
thermocouple  beads  were  soldered  into  the  grooves  with  a 
standard  lead- tin  low- temperature  solder.  The  leads  were 
secured  to  the  condenser  by  banding  with  twisted  wire  and 
were  passed  through  the  small  condenser  flange  to  the  mercury 
slip-ring  unit  on  t.he  test  stand.  Two  additional  thermo- 
couples sheathed  in  0.0625  inch  stainless  steel  tubing 
penetrated  the  main  flange.  The  sheathes  protruded  2.0  inches 
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into  the  evaporator  to  measure  the  vapor  temperature  and 
were  soldered  closed  at  the  tip  to  maintain  the  system 
vacuum.  The  inlet  and  outlet  temperatures  of  the  cooling 
water  were  monitored  by  Kapton  insulated  copper-constantan 
thermocouples.  All  thermocouples  were  read  directly  on  a 
Newport  Model  267  digital  pyrometer.  Calibration  was 
performed  in  a Rosemount  constant  temperature  bath  using  a 
platinum  resistance  thermometer  as  a standard,  A graph  of 
thermocouple  error  as  a function  of  thermocouple  reading 
was  plotted  to  be  used  to  correct  each  experimental  reading. 
Cooling  water  flow  rate  was  measured  with  a standard  rotameter. 
Accurate  regulation  of  the  flow  rate  was  accomplished  by 
adjusting  the  pressure  regulator  installed  in  the  cooling 
line . 

B.  EXPERI.MENTAL  PROCEDURES 

1 . Preparation  of  Condenser  Wall 

Tucker  [20]  established  the  procedure  for  prepara- 
tion of  the  copper  condenser  for  obtaining  film  condensation 
with  water.  This  procedure  was  followed  with  the  modifica- 
tions listed  below: 

a)  The  heat  pipe  was  tilted  so  that  the  evaporator 
was  depressed  slightly. 

b)  The  standard  cleaning  steps  were  followed  by 
scrubbing  with  a bristle  brush.  By  depressing 
the  evaporator  end,  cleaning  fluids  were  prevented 
from  running  back  into  the  heater  assembly. 
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c)  Complete  wetting  of  the  condenser  wall  was 
observed  after  cleaning. 

2 . Filling  and  Venting  Procedure 

Instead  of  using  the  filling  procedure  of  previous 
studies,  a simpler  procedure  was  used.  The  condenser  was 
depressed,  to  nearly  the  vertical,  with  the  evaporator  end 
raised  up.  The  heat  pipe  was  then  filled  with  about  350  ml 
of  working  fluid.  The  0-ring  and  glass  viewports  were 
installed.  The  bolts  of  the  retainer  ring  were  torqued 
to  thirty  inch-pounds  in  ten  inch-pound  steps  to  avoid 
cracking  the  glass  window.  This  filling  procedure  was 
followed  by  a venting  procedure,  used  by  Tucker  [20],  with 
the  modifications  listed  below: 

a)  Tilt  the  heat  pipe  with  the  evaporator  end 
down,  until  the  working  fluid  covers  up  the 
viewport . 

b)  Rotate  the  heat  pipe  by  hand  so  that  the  vent 
screw  is  on  top. 

c)  Slowly  apply  the  power  to  the  heater  to  about 
0.25  KW  until  the  saturation  temperature  in 

the  heat  pipe  is  at  least  220  degrees  Fahrenheit 
for  the  distilled  water. 

d)  Back  out  the  vent  screw  allowing  some  working 
fluid  vapor  to  escape.  Control  the  flow  of  the 
vapor  with  the  vent  screw,  being  careful  to  keep 
the  pressure  in  the  heat  pipe  above  atmospheric 
pressure . 
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e)  Vent  for  a minimum  of  ten  minutes  for  distilled 
water,  to  allow  escaping  steam  to  drive  off 
air  trapped  within  the  system,  leaving  about 
250  ml  of  working  fluid  left  over. 

f)  Turn  off  the  electrical  power  and  reseat  the 
vent  screw. 

g)  Return  the  heat  pipe  to  the  horizontal  position, 
and  turn  on  the  cooling  water  to  the  condenser. 

3 . Run  Procedure 

With  filling  and  venting  completed,  the  system  was 
ready  to  be  run.  The  standard  run  procedure  was  as  follows: 

a)  Open  cooling  water  supply  to  the  condenser  and 
bearing  cooling  water  lines. 

b)  Start  the  cooling  water  supply  pump,  and  adjust 
the  rotameter  to  the  desired  flow  rate.  In  this 
experiment,  only  the  50  percent  flow  rate  was  used 

c)  Oil  the  evaporator  bearing. 

d)  Energize  and  start  the  drive  motor.  Increase  the 
speed  to  1100  P.PM  to  insure  an  annulus  has  formed 
in  the  evaporator. 

e)  Set  rotational  speed  to  the  desired  value. 

f)  Let  the  system  reach  steady  state;  take  zero 
power  thermocouple,  cooling  water  flow  rate  and 
tachometer  readings. 

g)  Set  power  to  heater  to  desired  level  and  allow 
the  system  to  come  to  steady  state.  Take 
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tachometer,  cooling  water  flow  rate,  thermo- 
couple, voltage  and  ampere  readings, 

h)  Repeat  step  g)  for  each  power  level  desired. 

i)  Repeat  step  f) , 

j)  Repeat  steps  e)  through  i)  for  each  rotational 
speed  desired. 

4 . Data  Reduction 

Data  were  reduced  by  correcting  the  recorded  thermo- 
couple readings  to  obtain  a true  reading  in  degrees  Celsius. 
Correction  was  accomplished  using  thermocouple  calibration 
curves.  The  Celsius  readings  were  then  converted  to  degrees 
Fahrenheit.  Percentage  flow  rate  was  converted  to  pounds 
per  hour  using  the  rotameter  calibration  curve.  Heat  trans- 
fer for  each  run  was  calculated  as  the  product  of  the  cooling 
water  mass  flow  rate,  the  difference  in  temperature  between 
cooling  water  outlet  and  inlet  and  the  specific  heat  of 
the  cooling  water  using  the  basic  equation: 

Q = m C ,AT  (40) 

P 

where 

m = cooling  water  mass  flow  rate,  Ibm/hr 

C = specific  heat  of  the  cooling  water,  taken 

^ as  one  BTU/lbm-deg  F 

AT  = the  temperature  difference  between  the 
cooling  water  outlet  and  inlet,  deg  F. 
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To  correct  for  heat  generated  by  friction  in  the 
Garlock  seals  and  by  viscous  dissipation  of  the  cooling 
water,  the  average  zero  power  heat  transfer  rate  was  sub- 
tracted from  the  measured  heat  transfer  rate  in  order  to 
give  a corrected  heat  transfer  rate. 
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IV.  DISCUSSION  OF  EXPERIMENTAL  RESULTS 


A.  EXPERIMENTAL  RESULTS 

The  various  experimental  runs  are  shown  in  Table  11. 

All  the  heat  pipes  in  this  thesis  were  run  at  700,  1400  and 
2800  RPM,  except  the  1.00  inch  diameter  cylindrical  condenser 
which  was  run  at  1800  and  2800  RPM  for  the  fourth  run  and 
1400  and  2800  RPM  for  the  fifth  run  due  to  erratic  readings 
of  most  of  the  thermocouples  at  the  low  rotation  speed. 

It  was  concluded  that  the  problem  e.xisted  in  the  slip-ring 
unit  and  was  due  to  dirty  mercury  in  the  slip-ring.  So 
after  the  fifth  run,  the  slip-ring  unit  was  cleaned  by 
taking  the  mercury  and  using  methanol  to  clean  the  inside. 
Nitrogen  was  blown  into  the  unit  to  dry  it  out.  New  mercury 
was  added  prior  to  the  sixth  run.  The  plots  of  heat  transfer 
rate  versus  saturation  temperature  are  referred  to  as  the 
performance  curves. 

Figures  17,  18,  and  19  show  the  results  of  different 
condensers.  In  all  cases,  the  heat  transfer  rate  increased 
with  increasing  RPM.  In  the  cases  of  the  rotating  cylindri- 
cal condenser,  the  condensate  flow  is  induced  not  by  centri- 
fugal force  but  by  a hydrostatic  pressure  gradient  established 
in  the  condensate  as  a result  of  its  variable  film  thickness. 
By  increasing  the  diameter  of  the  cylinder  from  1.00  inch 
to  1.46  inches,  the  heat  transfer  capacity  increases.  For 
the  0.5  degree  truncated  cone,  the  heat  transfer  is  higher 
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TABLE  11 


Summary  of  Experimental  Runs 


Run 

Geometry 

RPM 

Test  Fluid 

No. 

1 

Cylindrical  1.46 

700 

Water 

inch  diameter 

1400 

2800 

2 

Cylindrical  1.46 

700 

Freon  113 

inch  diameter 

1400 

2800 

3 

Cylindrical  1.46 

700 

Ethanol 

inch  diameter 

1400 

2800 

4 

Cylindrical  1.00 

1800 

Water 

inch  diameter 

2800 

5 

Cylindrical  1.00 

1400 

Water 

inch  diameter 

2800 

6 

0 . 5 degree  trun- 

700 

Water 

cated  cone 

1400 

2300 

7 

0 . 5 degree  trun- 

700 

Water 

cated  cone 

1400 

2800 
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Figure  17. 


Heat  Transfer  Rate  Versus  Saturation  Temperature 
for  Cylindrical  1.00  Inch  Diameter  Condenser 
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90  100  110  120  130  140  150  160  170  180  190  200  210  220  230 

Saturation  Temperature,  deg  F 


Figure  18.  Heat  Transfer  Rate  Versus  Saturation  Temperature 
for  Cylindrical  1.46  Inch  Diameter  Condenser 
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than  the  cylindrical  condenser.  This  improvement  is  due 
to  the  centrifugal  force  which  increases  the  axial  force 
component  on  the  condensate  and  thereby  increases  the  con- 
densate flow  rate.  This  increase  in  working  fluid  circula- 
tion simply  means  that  more  vapor  can  be  condensed  per  unit 
time.  The  increase  in  centrifugal  force  also  flattens  the 
condensate  film  and  thereby  reduces  the  overall  thermal 
resistance  in  the  condenser. 

The  wall  temperature  profiles  as  shown  in  Figures  20, 

21  and  22  were  taken  from  the  1.46  inch  diameter  cylin- 
drical condenser  at  different  RPM.  Figure  23  was  taken 
from  the  0.5  degree  truncated  cone  at  700  RPM.  Distance 
zero  corresponds  to  the  condenser  end  nearest  the  evaporator. 
The  temperatures  would  be  expected  to  be  nearly  constant 
if  the  heat  pipe  were  truly  isothermal.  For  the  cylindri- 
cal condenser,  the  highest  temperature  is  at  the  evaporator 
end,  and  the  lowest  temperature  is  at  about  0.25  inch  from 
the  evaporator  end.  For  the  truncated  cone,  the  tempera- 
ture profile  decreases  almost  linearly,  and  the  highest 
temperature  is  at  the  evaporator  end.  This  increase  in 
temperature  at  the  evaporator  is  probably  due  to  conduction 
effects  from  the  heater  and  the  main  bearing. 

B.  ANALYSIS  OF  THE  RESULTS 

To  compare  the  theoretical  and  experimental  results, 
a log-log  plot  of  the  mean  Nusselt  number  versus  the  Sherwood 
number  was  made.  To  make  this  comparison  of  theoretical 


Figure  20.  Condenser  Outside  Surface  Temperature  Profile  at  700  RPM  for 
1.46  Inch  Diameter  Cylindrical  Condenser 


Distance  along  Condenser,  inch 

Condenser  Outside  Surface  Temperature  at  1400  RPM  for 
1.46  Inch  Diameter  Cylindrical  Condenser 


150 
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nser  Outside  Surface  Temperature  Profile  at  2800  RPM 
,46  Inch  Diameter  Cylindrical  Condenser 


Distance  along  Condenser,  inch 

Figure  23.  Condenser  Outside  Surface  Temperature  Profile  at  700  RPM 
for  0.5  Degree  Truncated  Cone 


versus  experimental  results,  the  following  simplifications 
were  made: 

1'  The  thermal  conductivity  of  the  wall  was  constant 
along  the  condenser, 

2)  the  inside  wall  temperature  was  found  by  using  a 
direct  application  of  Fourier's  law, 

3)  all  fluid  properties  were  evaluated  at  the 
saturation  temperature, 

4)  heat  flux  through  the  condenser  was  considered 
uniform, 

5)  for  the  truncated  cone,  the  surface  area  used  the 
average  radius  between  the  entrance  and  the  exit 
radius,  and 

6)  for  h = Q/AAT,  AT  was  taken  as  the  difference 

m 

between  the  saturation  temperature  and  the  condenser 
inner  wall  temperature. 

The  mean  Nusselt  number  for  a given  heat  transfer  rate 
is  given  by: 


Nu  = h L/k 
m m 


(Q/A  AT)L 
k 


(41) 


where 


Nu  = mean  Nusselt  number 
m 

2 

h = mean  heat  transfer  coefficient,  BTU/hr-ft  -deg  F 
m 

2 

A = the  average  surface  area,  ft 
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k 


thermal  conductivity  of  the  fluid, 
BTU/hr-ft-deg  F 


L = condenser  length,  ft 

The  mean  Nusselt  number  was  determined  from  the  data, 
and  this  value  was  plotted  against  the  corresponding  Sherwood 
number.  Figure  24  shows  Nu  versus  Sherwood  number  for 
the  cylindrical  condensers.  Fqir  water  with  the  1.46  inch 
diameter  cylindrical  condenser,  the  experimental  data 
agree  with  equation  (11).  The  data  using  the  1.00  inch 
diameter  cylindrical  condenser,  however,  show  that  the  Nu 

m 

increases  with  decreasing  Sh.  This  trend  is  the  opposite 
to  what  was  expected.  The  reason  for  this  trend  may  be 
due  to:  (1)  ripples  in  the  condensate  film  due  to  higher 

vapor  velocities  or  (2)  a high  uncertainty  in  the  low 
power  runs.  For  Freon  113  and  ehtanol,  with  the  1.46  inch 
diameter  cylindrical  condenser,  the  experimental  results 
are  higher  than  but  are  parallel  to  equation  (11) . The 
reason  for  this  discrepancy  is  not  clear  at  this  time.  It 
may  be  due  to  differences  in  vapor  densities  between  water. 
Freon  113  and  ethanol,  leading  to  differences  in  non- 
condensible gas  migration,  or  to  differences  in  vapor 
velocity  effects  within  the  theoretical  model. 

Figure  25  shows  the  mean  Nusselt  number  versus  the 
Sherwood  number  of  the  truncated  cone  to  find  out  the  best 
correlation  between  the  theoretical  and  experimental  results. 
The  results  of  Dhir  and  Lienhard  [8],  and  Roetzel  [14]  were 
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Figure  24.  Nusselt  Number  Versus  Sherwood  Number  for  Cylindrical  Condensers 


& Lienhard  [8] 


Figure  25.  Nusselt  Number  Versus  Sherwood  Number  for  Truncated  Cone  Condensers 


chosen  to  be  compared  with  the  experimental  results.  For 
the  one  degree  truncated  cone,  using  Tucker's  data  [20] 
and  Wagenseil's  data  [19],  the  experimental  results  for 
water  fall  between  both  analyses,  whereas  for  ethanol  and 
Freon  113,  the  experimental  results  are  higher  than  both 
analyses.  This  may  be  due  again  to  differences  in  vapor 
densities  between  water.  Freon  113  and  ethanol,  leading  to 
differences  in  non-condensible  gas  migration,  or  to  differ- 
ences in  vapor  velocity  effects  within  the  theoretical 
model.  For  the  0.5  degree  truncated  cone  with  water  as  the 
working  fluid,  the  experimental  results  fall  lower  than  both 
analyses,  but  fall  higher  than  the  analysis  of  Nimmo  and 
Leppert  [7]  for  a straight  cylinder. 

Figure  26  shows  Nu^/Nu_^  plotted  versus  the  parameters 
of  Roetzel's  analysis.  By  using  all  e.xperimental  results 
taken  to  data  at  the  Naval  Postgraduate  School,  it  appears 
that  good  agreement  exists  between  Roetzel's  analysis  and 
experimental  results.  The  differences  of  experimental 
results  and  Roetzel's  for  the  cylindrical  condensers  are 
less  than  15  percent.  For  the  truncated  cone  condensers, 
the  differences  are  less  than  20  percent  for  water.  For  the 
truncated  cone  with  ethanol  and  Freon  113,  experimental 
results  are  higher  than  Roetzel's  analysis  by  about  35 
percent.  However,  the  agreement  is  reasonable  and  the 
analysis  of  Roetzel  can  be  used  to  predict  rotating  heat 
pipe  performance  to  within  i 20  percent  for  water.  For  Freon 
113  and  ethanol  the  data  lies  above  the  theory  by  +35  percent. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

For  an  internally  finned  condenser,  the  analysis  of 
the  condenser  heat  transfer  problem  using  a two-dimensional 
wall  conduction  model  and  a parabolic  temperature  distribution 
along  the  fin  can  be  approximated  by  the  two-element  Finite 
Element  solution  with  twenty-five  axial  increments. 

For  a smooth  condenser; 

1.  The  theoretical  work  of  Leppert  and  Nimmo  [6,7]  can 
be  used  to  predict  satisfactorily  cylindrical  heat 
pipe  performance. 

2.  The  theoretical  equation  which  best  predicts  overall 
heat  pipe  performance  is  due  to  Roetzel  [14] : 

Nu^  = Nu  [1  + 1.924  tan  $ . 

00 

The  above  equation  agrees  with  the  experimental  data 
of  water  to  within  t 20  percent.  The  experimental 
data  of  Freon  113  and  ethanol  lie  approximately  35 
percent  above  the  theory. 

B.  RECOMMENDATIONS 

Continue  to  study  the  internally  finned,  rotating  heat 
pipe  by: 
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1.  Re-examining  the  Finite  Element  analysis  using  the 
mirror  geometry  presented  in  Figure  7. 

2.  Analyzing  film  condensation  on  a variety  of  fin 
shapes,  and  applying  the  solution  of  the  two- 
dimensional  model  to  determine  the  shape  permitting 
optimum  heat  transfer  rate, 

3.  Manufacturing  an  internally  finned  condenser  to 
obtain  experimental  data  to  compare  with  the 
theoretical  analysis. 

Continue  to  study  the  smooth,  rotating  heat  pipe  by: 

1.  Using  a condenser  section  with  a half  cone  angle 
greater  than  1 degree  to  obtain  additional  experi- 
mental data  to  compare  to  Roetzel's  analysis. 

2.  Examining  reasons  why  the  Freon  113  and  ethanol 
data  are  falling  above  the  theory. 


